Casimir effect due to fluctuations of geometry in a quantum antiferromagnet 
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We show that fluctuations of geometrical connectivity can give rise to a Casimir type attractive 
or repulsive force in a two dimensional Heisenberg antiferromagnet. The type of fluctuation that 
we consider, called a phason flip, is well known in quasicrystals, but less so in periodic structures. 
As the classical ground state energy of the antiferromagnet is unaffected by this type of fluctuation, 
energy changes are purely of quantum origin. We calculate the effective interaction between two 
parallel domain walls, defining a slab of thickness d, in such an antiferromagnet within linear spin 
wave theory. The interaction is anisotropic, and for a particular orientation of the slab we find that 
it decays as 1/d, thus, more slowly than the electromagnetic Casimir effect in the same geometry. 
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FIG. 1. The staggered dice lattice. The sites are colored 
red or blue according to whether they belong in the A or B 
sublattice. 



The Casimir effect [l| refers to the net force between 
two electrically neutral objects that arises due to vac- 
uum fluctuations of the electromagnetic field when they 
are placed suffiently close together. Initially considered 
as a somewhat mysterious apparition of "something out 
of nothing" , the phenomenon described by Casimir is far 
more widespread, and arises quite generally due to fluc- 
tuations in media where long ranged correlations [2[ are 
present. Thermal fluctuations in superfluids 3] and su- 
perconductors [J] give rise to long ranged forces. Casimir 
forces have been shown to exist in nematic liquid crys- 
tals or even granular systems @. In this paper we 
consider the interaction energy in a thin slab of the spin- 
i antiferromagnetic Heisenberg model on the staggered 
dice lattice illustrated in FigQ] The interfaces corre- 
spond to a row of phason flips of the structure, and hence 
to boundary conditions of a novel type compared to the 
usual Dirichlet or Neumann boundary conditions. 

In our model, periodic boundary conditions are as- 
sumed so that there are no external boundaries. The 
parallel slab geometry in this antiferromagnet is instead 



determined by geometrical defects termed phason flips, in 
which the bonds of a single spin flip between two equiva- 
lent local configurations. We assume that the transitions 
between the two configurations can occur easily and that 
the magnitude of the nearest neighbor Heisenberg spin- 
spin couplings, J, is unaffected by the change. Such a 
situation could arise in structures in which there are two 
equivalent chemical pathways to generate antiferromag- 
netic Heisenberg superexchange, as we will discuss at the 
end of this paper. This type of disorder is reminiscent 
of the "two level systems" that were proposed for amor- 
phous systems Q, in order to explain the anomalously 
high scattering rate of phonons seen at low temperatures. 
While phason flip disorder is well-known in the field of 
quasicrystals [8j, its effects have only recently been con- 
sidered in a periodic structure, namely the staggered dice 
antiferromagnet [9]. Compared to the disordered quan- 
tum spin models studied in the literature such as bond- 
disordered models or site-dilution models ljl|), pha- 
son flip disorder leads to some novel properties. In con- 
trast to bond-fluctuation or site-dilution, a phason flip 
does not modify the classical ground state energy of the 
antiferromagnet. It does however modify the quantum 
corrections to the ground state energy and the spectrum 
of excitations around the ground state. As a result, an- 
tiferromagnetism was found to be enhanced rather than 
diminished as in other types of disordered antiferromag- 
nets . This phenomenon can be considered as a variant 
of the order-by-disorder effect originally found by Villain 
for a frustrated [HI spin model [Hi ]. 
The Heisenberg model considered is 
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with J > 0, where Si are the spin operators for the sites 
i of the staggered dice lattice and the sum runs over all 
pairs of spins i,j joined by an edge of length a. On 
the dice lattice, one can distinguish hexagonal "cages" 
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formed by six spins, of which three are coupled to the spin 
at the cage center. This can occur in one of two ways. 
The staggered dice lattice is formed by assembling such 
hexagonal cages such that the bonds alternate between 
the "up" and "down" configurations. The lattice vectors 
of the SDL shown in Figfj] are given by r\ = \fiax and 
r-2 = 'Say, and the symmetry of the structure is rectan- 
gular. The SDL structure is unfrustrated, as can be seen 
in FigJTJ where spins are colored red or blue, according 
to the sub-lattice and one expects rotational symmetry 
to be broken in the ground state, which is Neel ordered, 
with spins on sublattice A aligned along, say, the +z- 
direction, and spins on sublattice B aligned along — z. 
We recall that in the ordinary dice (also known as T3) 
lattice, the central sites are linked to three neighbors in 
a uniform way instead of this alternating st agg ered fash- 



ion, and that the result is a ferrimagnet [171 ]. In spin 
wave theory, one represents the spin operators in terms 
of Holstein-Primakoff boson operators a,(aj) for sites i on 
the A-sublattice and bi(b]) for sites i on the B-sublattice 
(ljj . The linear spin wave Hamiltonian for the SDL lat- 
tice of N sites obtained from EqfT] after neglecting higher 
order terms is 
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One can introduce the Fourier transformed set of opera 
tors 
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where the sum runs over the N c = N/6 unit cells sit- 
uated at positions R n and /i = 1,2,3 labels the three 
sites belonging to each sublattice within a unit cell. The 
Hamiltonian now reads 

=JSJ2J2 z » V ( a U a ^ + b^) + 

k 

(7mm' (k)a^k Vfe + h - c -) ( 5 ) 

where — 4, 5 and 3 are the coordination numbers for 
/i = 1, 2, 3 and j(k) is the 3x3 matrix 
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with — exp(ik.ffj,). We consider henceforth the solu- 
tions obtained for S = h, the case for which the quan- 
tum fluctuations are strongest. The energies of the spin 
excitations in the model are determined upon diagonal- 
ization by means of a Bogoliubov transformation. In the 




FIG. 2. Bond configurations before (left) and after (right) 
two phason flips. The six sites of the unit cell are numbered. 



limit that k tends to zero one finds the linear relation- 
ship for the lowest energy Goldstone modes tuj(k) = v s k 

(for k <C a -1 ) where v s = G^J^J is the spin wave veloc- 
ity. The ground state energy per site was calculated in 
linear spin wave theory to be Eg = —0.6517(7) J in the 
thermodynamic limit, which compares well to the result 
Eq = —0.6639(1) J obtained from quantum Monte Carlo 
(QMC) simulations [§]. 

Fig [2] illustrates the effect of a pair of phason flips 
within a single unit of the dice lattice. The interior spins 
S3 and Sq are coupled with three cage spins in different 
ways in the left hand and the right-hand figures. As the 
total number of bonds is constant, the total classical en- 
ergy is unchanged by phason flips. As can be seen from 
Fig|21 the transition between the two configurations im- 
plies a change of the sublattice of the center spin. In 
order to preserve the antiferromagnetic condition on the 
total spin, Stot = 0, the number of flips on sublattices A 
and B respectively, are therefore taken to be equal in all 
the (finite) systems that we consider. Since phason flips 
do not introduce any frustration, Neel antiferromagnetic 
order is always preserved and linear spin wave theory can 
thus be used to study the system within the entire range 
of values of the number of flips, < Nf < N/3 where N 
is the total number of spins. Periodic boundary condi- 
tions are assumed along x and y directions. 

The energy cost per phason flip can be calculated in 
spin wave theory by a numerical diagonalization proce- 
dure. It was found to be approximately 0.06 J, far smaller 
than, for ex amp le, the energy of a spin vacancy which is 
about 0.6 J [19(. The energy of interaction as a function 
of the distance has a short distance component - of the 
order of one unit cell spacing - and an anisotropic power 
law decay at long distance. The sign of the interaction 
energy of two phasons depends on the sublattice: for pha- 
sons on opposite sublattices, the interaction is attractive, 
and for those on the same sublattice it is repulsive. The 
rapid decay of perturbations around a phason flip can be 
understood in terms of a continuum version of the dis- 
crete Hamiltonian of Eq[T]by the fact that a flip produces 
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FIG. 3. Illustration of domain wall configurations used in 
the study of Casimir effect, a) Slab aligned along fi direction 
(shaded light blue) b) Slab aligned along fj + r*2 direction. 
Sites on the interfaces belong in either the A-sublattice (red 
points) or B-sublattice (blue points). 



gradient terms of high (third and above) order. 

We now present results for the energy dependence as 
a function of thickness for a horizontal (i.e. having its 
orientation parallel to the r\ direction) slab of thickness 
d. Fig,a shows the configuration studied: the lower row of 
A-sublattice phasons is shown in red, while the upper row 
corresponds to B-sublattice phasons, shown in blue. The 
resulting system is analogous to a parallel plate capacitor: 
outside" the parallel lines one has the original dice lattice, 
while "inside" the plates one has the fully flipped version 
of the original lattice. For each value of d, we calculated 
the difference of ground state energy E^(d) — Eq where 
E is the ground state energy of the defect-free SDL. 
Dividing this by the length of the domain wall L X1 we 
obtain the energy per unit length, e. As Fig0] shows, the 
interaction between the lines is attractive. We fitted the 
energy change to the expression 
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where the second term of the term in brackets arises due 
to the periodic boundary conditions (L y is the periodic 
length in the y-direction). Fitting to the data using EqJT] 
for distances that are sufficiently large, we obtain the co- 
efficients An and Bn for each sample size. The energy 
versus distance plot for three different sample sizes are 
shown in FigJU along with the curves given by Am+Bm / d 
for N = 4800, 5888 and 6912. Both A and B coefficients 
follow a finite size scaling in 1/N, in accordance with the 
scaling expected for the energy per length of a ID anti- 
ferromagnet. Extrapolation to the infinite system yields 
the interface energy per length: e(d) = A + B/d, with the 
values A = 0.062 J/a and B = 0.001J. If one takes the 
exchange coupling J to have a value comparable to that 
of typical Heisenberg antiferromagnets, namely a frac- 
tion of an eV, and that the edge length is of the order of 
an Angstrom, one obtains an interaction energy of about 
lpj/m. 




FIG. 4. Change in ground state energy per unit length e(d) 
plotted as a function of the distance d between two rows of 
defects for three system sizes (N = 4800,5808,6912). The 
lines represent a fit to Eq(7] The inset shows in a log-log plot 
the scaling with 1/iV of the An and Bn coefficients. 



When the horizontal slab calculations are repeated but 
for domain walls of phason flips on the same sublattice, 
the resulting force is repulsive, with the same d depen- 
dence. This behavior is consistent with the earlier stud- 
ies of the interaction between two phason flips: repul- 
sive for same sublattice, attractive for different sublat- 
tices Q . This mechanism for repulsive Casimir force can 
be contrasted with the system described in [2(| wherein a 
topological insulator is used as a dielectric medium, and 
in which the sign of the interaction is varied by tuning 
the magnetoelectric polarizability of the medium. The 
Casimir effect is strongest for the horizontal configura- 
tion, where the domain walls are rows of closely spaced 
same-sublattice flips. All other orientations of the do- 
main walls lead to faster decay. In particular, a domain 
wall along the r\ +r2 direction for example (Fig.), which 
alternates flips on the A- and B-sublattices, leads to a 
short range contact-type interaction. 

In sum, the interaction energy E{d) between domain 
walls which are parallel to x depends on the separation d 
as —L x /d. This dependence on the inverse of the thick- 
ness can be guessed from a dimensional analysis. On 
the one hand, the interaction energy is a dimensionless 
quantity when expressed in units of J (the only energy 
scale in the problem), and on the other hand, this en- 
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ergy is expected to be proportional to the length of the 
interface, which the interaction energy is proportional 
to L x /d. This distance dependence of the energy can 
be compared with that of the electromagnetic field in 
a similar geometry. In that dimensional analy- 

sis predicts that the interaction energy due to quantum 
fluctuations, which is proportional to he, should drop 
off as L x /d 2 . Our model also differs in several essen- 



tial respects from the one considered in [2l| between two 
rows of holes in a square lattice antiferromagnet, where 
a Casimir effect is observed. For that problem, the in- 
teraction, always attractive was found to decay as 1/d 2 . 
In that case, each row of holes is characterized by the 
microscopic length scale a, and the region between is the 
unperturbed antiferromagnet. In contrast, the slab ge- 
ometry that we consider separates two regions which are 
distinct, although related by phason flips. 

A back-of-the-envelope calculation for this inverse dis- 
tance dependence of E can be given. We note that at long 
distance, the two phason interaction energy decays as r~ 3 
@ • The d dependence of the energy of the slab of phasons 
can be estimated by computing the energy of interaction 
of a single A-sublattice phason on the edge due to its in- 
teractions with successive rows of B-sublattice phasons, 

which is approximately Jq dy Jq"^ 2 dx(x 2 + y 2 )~ 3 ^ 2 ~ 
1/d for L x ^> d ^> a. Added to this is a contribution of 
the opposite sign, due to the interactions of same sublat- 
tice phasons, but the distances are different and therefore 
the contributions will not completely cancel. This rough 
calculation for the energy per unit length makes use of 
the simplified asymptotic form of the 2-phason interac- 
tion, and does not take into account the detailed angular 
dependence of the anisotropic decay. 

We briefly address the question of possible experimen- 
tal realizations of the ordinary dice and staggered dice 
lattices. On the one hand, the ordinary dice lattice can 
been realized using optical trapping of cold atoms [22[ , as 
a realization of tight binding models with so-called Dirac- 
Weyl fermions. It would be interesting to investigate pos- 
sibilities for using such systems to realize the Heisenberg 
model. On the other hand, in conventional condensed 
matter systems, the staggered dice lattice which is a rel- 
atively complex structure, could perhaps be realized in 
layered compounds. As an example of how this may be 
done, we can consider how the Heisenberg model on the 
two dimensional Kagome lattice (dual of the dice lattice) 
has been experimentally realized in a number of com- 
pounds. One of these, called volborthite, has spins asso- 
ciated with Cu atoms, which interact via super-exchange 
processes. When each of these magnetically active lay- 
ers is viewed from above, the spins and the bonds are 
seen to project to the Kagome structure. In an analo- 
gous way, our model on the dice lattice could be realized 
by considering spins on vertices of cubes, and considering 
their projection on a plane perpendicular to one of the 
body diagonals that bisects the cube. The six sites of 





FIG. 5. (upper) Two examples of rows of cubes showing the 
vertices before projection. Red sites project onto hexagonal 
cage vertices, while the blue site can lie either above or below 
the 3-fold plane bisecting the cubes. The upper row shows a 
uniform choice of blue sites, and the lower row a staggered 
choice. (lower)The sites after projection on the three-fold 
plane. The upper row has the same decoration in all cells, 
while the lower row alternates between two configurations. 



the hexagonal cell of the dice lattice correspond to the 
six vertices lying closest to this plane (shown in red in 
Fig E] for a row of cubes). The two remaining apical ver- 
tices shown in blue, both project onto the central site 
in the cage but with different connectivities. If the one 
of these configurations is selected randomly, for example, 
one obtains the phason-disordered case. One could in 
principle impose a given configuration by suitably engi- 
neering the substrate on which the spin system is built, so 
as to slightly bias one or other of the two configurations. 
Transitions could then be induced by small changes of 
the distances of the apical atoms. 

In conclusion, we have computed the Casimir forces in 
a two dimensional antiferromagnet based on the honey- 
comb lattice. The fluctuations have a geometrical origin, 
involving changes of local topology, and may be realized 
in quasi-two-dimensional structures in which the nearest 
neighbor environments allow for local transitions between 
two different conformations with little or no energy cost. 
The interaction energy depends on the orientation of the 
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Casimir slab, and can be attractive or repulsive depend- 
ing on whether the walls belong to the same to different 
sublattices. The effect is strongest for a slab with hori- 
zontal orientation, the interaction energy falling off as the 
inverse of the slab thickness. The magnitude of the effect 
is estimated to be of the order of pJ/ m, for antiferromag- 
netic couplings J of the order of leV. This fluctuation- 
induced interaction could be measured for sheets of spins 
in which the phason flips are flips between isomeric con- 
figurations. 
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